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We analyze the physical properties of boson stars, which possess counterparts in flat space-time, 
Q-balls. Applying a stability analysis via catastrophe theory, we show that the families of rotating 
and non-rotating boson stars exhibit two stable regions, separated by an unstable region. Analogous 
to the case of white dwarfs and neutron stars, these two regions correspond to compact stars of lower 
and higher density. Moreover, the high density phase ends when the black hole limit is approached. 
Here another unstable phase is encountered, exhibiting the typical spiralling phenomenon close to 
the black hole limit. When the interaction terms in the scalar field potential become negligible, the 
properties of mini boson stars are recovered, which possess only a single stable phase. 
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I. INTRODUCTION 

After stars have consumed their nuclear fuel, they end their lives as compact astrophysical objects, and turn into 
white dwarfs, neutron stars (or their variants) or black holes, depending on their initial mass. When one investigates 
the equilibrium properties of compact stellar configurations, as discussed for instance in Shapiro and Teukolsky 
one finds two stable phases for such compact stellar objects. These stable phases correspond to the white dwarf 
phase and the neutron star phase, where equilibrium is achieved by the electron and the neutron degeneracy pressure, 
respectively. The two stable phases are separated by an intermediate unstable phase. Moreover, the stable neutron 
star phase is followed by an unstable phase, exhibiting a spiralling behavior (when the mass is considered as a function 
of the radius). As seen in Fig. 1, the stable neutron star phase ends and the unstable phase sets in, when the black 
hole limit is approached. 
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Figure 1: Schematic diagram of the mass versus the radius (left) and the mass versus the central density (right) of compact 
stars after Shapiro and Teukolsky 



Representing hypothetical astrophysical objects, boson stars have received considerable attention since being pro- 
posed by Feinblum and McKinley @], Kaup 0, and Bonazzola and Ruffini Boson stars are obtained when a 
massive complex scalar field is coupled to gravity. The conserved Noether current associated with the global U(l) 
symmetry is related to their particle number Q. The physical properties of boson stars depend crucially on the rest 
mass of the bosons and on the presence and type of self- interaction of the bosons (see e.g. the review articles 043)- 

When only a mass term is present but no self-interaction of the scalar field, the resulting family of mini-boson stars 
possesses one stable branch of solutions, which ends at the maximal mass configuration and is followed by an unstable 
set of solutions with spiralling behaviour. The maximal mass M max is on the order of the Planck mass squared, 
divided by the boson mass: Mpj/ms- When a repulsive quartic self- interaction is included the resulting family of 
boson stars retains a single stable branch, but the maximal mass M max is now on the order of yXM^/m^, where A is 
the self-coupling constant [lCj |. Thus as shown by Colpi, Shapiro, and Wasserman [l(| boson stars with much larger 
masses are obtained. 

Friedberg, Lee and Pang [Tljj considered boson stars with a self-interaction of the scalar field that allows for non- 
topological soliton solutions, also called Q-balls, even in the absence of gravity. They estimated that the maximal 
mass of such boson stars is on the order of Mpj/mg 
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Boson stars are obtained as stationary solutions of the coupled Einstein-scalar field equations, when the scalar 
field has a harmonic time dependence. The associated frequency w s is bounded from above by the scalar mass me, 
while its lower bound depends on the details of the model. Boson stars may rotate (l2T[l8l ] . Since their total angular 
momentum is quantized in terms of their particle number, J — nQ, boson stars may change their angular momentum 
only in discrete steps [l2l ] . 

The stability of boson stars has been addressed from various points of view. Whereas Lee and Pang performed 
a linear stability analysis of boson stars with respect to small oscillations, Kusmartsev, Mielke and Schunck [20, 21] 
applied catastrophe theory to extract the stable branches of families of boson stars. Catastrophe theory had been 
introduced by Rene Thorn in the 1960's [22j |. with further developments made by Zeeman (23|, Poston and Stewart 
[24], [25[, Arnol'd [2{||27|, and many others. Applications to solitons had been discussed by Kusmartsev [28j . 

Here we analyze the properties and in particular the stability of boson stars obtained with quartic and sextic self- 
interaction terms, as introduced by Friedberg, Lee and Pang in their non-topological soliton model [lll[3l[. Indeed, 
the existence of a flat space-time limit is of profound importance for the properties of these boson stars, as we will show 
below. With respect to stability we adapt the procedure of Tamaki and Sakai [2!| [3(| ■ They analyzed the stability of 
spherically symmetric Q-balls and boson stars via catastrophe theory under primarily mathematical points of view. 
Here we extend this analysis to rotating boson stars. 

In Section II we present the action, the equations of motion and the definition of the global charges for non-rotating 
and rotating boson stars. In Section III we analyze the physical properties of the families of boson stars, presenting 
first the results for the non-rotating case and then for the rotating case. We give our conclusions in Secion IV. The 
Appendix addresses briefly the construction of the solutions. 



II. MODEL 
A. Action 



We consider the action of a self-interacting complex scalar field $ coupled to Einstein gravity 

R 1 



S = 



-gd x , 



(1) 



16ttG 2 : ' 1 " 

where R is the curvature scalar, G is Newton's constant, the asterisk denotes complex conjugation, 

<j> = 

and U denotes the potential 

C/(|$|) = A|$| 2 (|$| 4 - a|$| 2 + b) = A (</> 6 - ac}) 4 + b<t> 2 ) , 

with |$| = (j>. The potential is chosen such that nontopological soliton solutions [llj . also referred to as Q-balls, exist 
in the absence of gravity. The potential has a minimum, U(0) — 0, at $ = and a second minimum at some finite 
value of |$|. The boson mass is given by me = yXb. 

Variation of the action with respect to the metric leads to the Einstein equations 



(2) 
(3) 



with k = 8ttG and stress-energy tensor T f 



9^ v Lm 



t dL_ 



Variation with respect to the scalar field leads to the matter field equation, 

dU 



□ 



am 



$ = 



(4) 

(5) 
(6) 

(7) 



where □ represents the covariant d'Alembert operator. Equations (j4j and jJJ represent the set of the coupled 
Einstein-Klein-Gordon equations. 

Also boson stars with two complex scalar fields have been considered, leading to interesting phenomena due to their 
interaction [32|, |33|, Hf| • The influence of a negative cosmological constant has been investigated in [34] . 
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B. Ansatz 



To obtain stationary axially symmetric solutions, we impose on the space-time the presence of two commuting 
Killing vector fields, £ and r], where 



(8) 



in a system of adapted coordinates {t, r, 9, ip}. In these coordinates the metric is independent of t and ip, and can be 
expressed in isotropic coordinates in the Lewis-Papapetrou form [36rl39| 



ds 2 = -fdt 2 + ~ 



h (dr 2 + r 2 d6 2 ) + 



r 2 sin 2 I 



\d<p dl 



(9) 



The four metric functions /, I, h and u are functions of the variables r and 9 only. 

The symmetry axis of the spacetime, where r\ = 0, corresponds to the z-axis. The elementary flatness condition 



x^x^ 

4X 



= 1 



X = 77% 



then imposes on the symmetry axis the condition [3a ] 

h\e=Q — h\ s=7r = 1 . 
For the scalar field <3? we adopt the stationary ansatz [l2[ 

<f>(t,r,6,<p) = <f>(r,0) e iu " t+in{p 
where <p(r, 9) is a real function, and uj s and n are real constants. Single- valuedness of the scalar field requires 

$(<p) = $(2tt + tp) , 



(10) 



(11) 



(12) 



(13) 



thus the constant n must be an integer, i.e., n = 0, ±1, ±2, .... We refer to n as rotational quantum number, since 
for n ^ axially symmetric rotating boson stars arise, whereas for n — spherically symmetric non-rotating boson 
stars are obtained. 

Solutions with positive and negative parity satisfy, respectively, 



(r,7r-0) = <j){r,9) 
(r,7T-0) = -<j>(r,9) 



(14) 
(15) 



C. Global charges 

The mass M and the angular momentum J of stationary asymptotically flat space-times can be obtained from their 
respective Komar expressions (40| . 

M =^q I R^CdV , (16) 

and 



J =-^G^ R ^ n ^dV . (17) 

Here S denotes an asymptotically flat spacelike hypersurface, n M is normal to E with n^n^ = — 1, dV is the natural 
volume element on S, £ denotes an asymptotically timelike Killing vector field and r] an asymptotically spacelike 
Killing vector field [4(|. The mass M and the angular momentum J can be read off directly from the asymptotic 
expansion of the metric functions / and u>, respectively, [38| 

M = lim r 2 a r / , J = -p^ lim r 2 u> . (18) 

A conserved charge Q is associated with the complex scalar field $, since the Lagrange density is invariant under 
the global phase transformation 

$ $e M (19) 
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leading to the conserved current 

f = -i - $9^$*) , j% = 0. 

The conserved scalar charge Q is obtained from the time-component of the current, 

f \g\ 1/2 drdOdip 



(20) 



(21) 



As first derived by Schunck and M ielke \12\ , one obtains a quantization relation for the angular momentum in terms 
of the charge, 



J = nQ . 

Thus a spherically symmetric boson star has angular momentum J = 0, because n = 0. 



(22) 



D. Units 

We choose for the potential U(<p), Eq. ©, the following set of fixed parameters [13, [HUHl 

A = 1, a = 2, b= 1.1 = (j% . 
The equations then depend only on the dimensionless coupling constant k 



k = 8ttG 



Mo 



(23) 



(24) 



Since k consists of a product of Newton's constant and the square of the boson mass, we may interpret a change of its 
numerical value in two ways: either as a change of the gravitational coupling for a fixed boson mass, or as a change 
of the boson mass for a fixed value of the gravitational coupling. 

To obtain the respective dimensionful values of the physical properties, M phys , Q phyS j J phys , we have to scale the 
numerically calculated values for the mass M num , the charge Q num and the angular momentum J num appropriately. 
Therefore we introduce the parameter go 



qo = 



where mp\ is the Planck mass and find 



The length scale is set by 



gphys 

M phys 
jphys 



Mo 



\ TUB J 87T 



^q Q aum , 
IJ,oKqomBM r 
nhQ nunl . 



Mo 



x 0.2 x 10" 15 m 



tib [GeV] 

Two examples for the scales of the physical properties are exhibited in Table 1. 



(25) 



(26) 



m B 


90 


qom B 


a 


lGeV/c 2 


f« 10 36 


w 10 36 GeV/c 2 


~ 0.2/iofm 


10- 19 GeV/c 2 


w 4 x 10 74 


«M sun 


w 0.2^ 10 4 m 



Table I: Scale factors for the physical properties 
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III. BOSON STAR PROPERTIES 



A. Spherically symmetric boson stars 



1. Equilibrium space 

Let us first consider the families of fundamental spherically symmetric boson stars (n = 0) as obtained in [l7l [Lsj . 
(In fundamental boson stars the scalar field is a monotonically decreasing function, whereas in radially excited boson 
stars, it possesses nodes.) For a fixed value of the dimensionless coupling constant k, Eq. (f24| . a family of stationary 
solutions exists in the frequency range ojq(k) < uj s < y m p. The minimal frequency u)q(k) depends on k and increases 



1 71 . The maximal frequency w max is always given by the boson 



with k, tending to finite limits as n — > and K — > oo 
mass. 

All families of solutions together form the equilibrium space, which we denote by M. = {w s , k, Q}. To illustrate the 
equilibrium space we have exhibited several families of solutions in Fig.[2j For each family (with fixed n) the frequency 
u> s is shown versus the particle number Q. The coupling constant k assumes values in the range 0.0002 < k < 1 in 
the figure. The equilibrium space may then be pictured as the surface obtained when k varies continuously from zero 
to infinity. 
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Figure 2: Equilibrium space M = {io 3 , Q, k} for fundamental boson stars in the range 0.0002 < k < 1. 

The families of solutions have two important features in common. They all start from particle number Q = at 
the upper limit of the frequency, and they all form spirals towards their lower frequency limit. Inbetween, however, 
a qualitative change of the curves is observed, as the value of k increases. The salient maximum value of Q in the 
small frequency range as well as the finite (relative) minimum value of Q, which are both present for small values of 
the coupling constant k, become less pronounced with increasing k, until they merge and disappear altogether. 

We may also consider the set of equilibrium solutions with respect to another set of parameters. Instead of the 
frequency oj s we may consider the finite value <f>o = 0(0) of the scalar field at the origin. (We recall, that in these 
fundamental boson star solutions the scalar field decreases monotically from the origin to zero at infinity.) Moreover, 
instead of the particle number Q we may consider the mass M of the solutions. 

We exhibit the same families of solutions as above for this choice of variables in Fig. [31 Here for each family of 
solutions (with fixed k) the value of the scalar field 4>o is shown versus the mass M. We note, however, that the range 
of 4>o has been truncated in the figure, being limited to the range < <fio < 2. Inside the spiral <fio certainly assumes 
larger values, and indeed increases monotonically as the mass exhibits damped oscillations (17| . The equilibrium space 
J\f = {(f>o, M, k} may then again be pictured as the surface obtained when k varies continuously from zero to infinity. 



2. Binding energy and cusp structure 



To get a better physical understanding of these boson stars that form the equilibrium space, let us next address their 
binding energy B = m^Q — M. We exhibit the binding energy B in Fig.|4]for two families of solutions, corresponding 
to typical examples in the lower k range, namely k = 0.01 and k = 0.02. While at first glance all solutions seem to be 
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k = 0.1000 




Figure 3: Equilibrium space M = {<j>o, M, k} for fundamental boson stars in the range 0.0002 < k < 1. 



bound, with their binding energy increasing almost linearly with their particle number, a closer look at the solutions 
in the high frequency range (uj s close to tub) reveals a bifurcation structure involving two cusps. 




200 



Figure 4: Binding energy B = M — tubQ for two families of solutions, corresponding to n = 0.01 and k = 0.02, respectively. 
The right figure zooms into the high frequency range. 

Clearly, a first branch of bound boson star solutions resides between the vacuum M = Q = 0, and a local maximum 
of the mass and charge, Ma, Qa, at the first cusp A. A second branch, where the mass and charge decrease 
monotonically, connects the first cusp A with the second cusp B, where the mass and charge have a local minimum, 
Mb , Qb- Most of this branch resides in the unbound region. The third branch emerges from the second cusp B and 
extends up to the third cusp C, where the mass and charge reach their global maximum, Mc, Qc- Along this third 
branch the mass and charge increase monotonically. The solutions are unbound only in the vicinity of the second 
cusp B, since the binding energy crosses the zero already at the point u, located closeby. At the third cusp C also 
the binding energy is maximal. 

We exhibit in Fig. [5] the physical characteristics of these 3 cusps, A, B and C. As k — > 0, the values of the mass 
Ma and the charge Qa at the cusp A increase, tending to infinity with the power k -1 / 2 , as illustrated in Fig. [5] by 
the curves MA.nm and QA.iim- The values of the mass Mb and the charge Qb at the cusp B show a very different 
behavior. They tend to constant values in the limit k — > 0, corresponding to the values of the mass and the charge of 
the unique minimum of the Q-ball solutions of flat space-time. This is also illustrated in Fig. [S] At the cusp C the 
values of the mass Mc and the charge Qc increase again without bound as k — > 0. As extrapolated previously [l7j . 
the limiting behaviour is of the form Mc ~ k~ 3 / 2 , Qc ~ k~ 3 / 2 , which is also shown in the figure. 

Turning to the larger values of K, we note that the two extrema B and C, and thus the two cusps B and C, merge 
and disappear at a critical value n CT . Thus beyond K cr , from the three cusps A, B and C only the cusp A is left. The 
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Figure 5: Coupling constant k versus the mass M (left) and versus the particle number Q (right) for the cusps A, B and C, 
together with functions approximating their behavior in the limit k — > 0. Note, that for the larger values of k the cusps B and 
C axe no longer present. 



large value of k means from a physical point of view, that the higher order terms in the potential become relatively 
less important, and the solutions tend to the well-known mini-boson star solutions. Indeed, for k — > oo the mini-boson 
star solutions are recovered after a rescaling (l7| . 

To complete the discussion of the branches of the families of boson star solutions, we still need to address their 
behavior beyond the cusp C (or for large k beyond the cusp A). As seen in Fig. [2] the branches form a spiral when 
considered in terms of the frequency w s and the charge Q (or the mass M). Alternatively, when considered in terms 
of the value of the scalar field at the origin and the mass M (or the charge Q) as in Fig. [31 a damped oscillation is 
seen. Thus beyond C (or for large k beyond A) the mass and the charge reach a minimum D, then another maximum 
E, another minimum F, etc., converging towards limiting values M\i m , Qn m [171 ]. When the mass is considered as a 
function of the charge, finally, this behavior translates into an intricate cusp structure [Tl| . 



3. Size and black hole limit 



There is no unique definition for the radius of a boson star, and many proposals have been discussed in the literature 
(see e.g. @). Here we have chosen the definitions 

J f\g\ 1/2 rdr 2 J f \g\ 1/2 r 2 dr 

1 — 372 ' 2 = T/2 \ ' 

Jf Iff I dr ' J f \g\ dr 

for the spherically symmetric boson stars, where the radial coordinate is not the isotropic coordinate introduced in 
Section II, but a Schwarzschild-like coordinate, which has an invariant circumferential meaning. Fig. [6] (left) shows, 
that both definitions give rather similar results for the size of the boson stars. From the radius and thus the size of 
the boson stars together with their mass we obtain an estimate of their density. 

Let us now discuss the relation between the mass and the size of the boson stars, focussing on the smaller values of 
k (away from the mini-boson star limit). As seen in Fig. El the mass increases along the first branch, while the radius 
decreases. This is the behavior known from compact objects such as white dwarfs and neutron stars. The first branch 
ends at the cusp A. Beyond A the mass decreases until the second cusp B is reached, while the radius continues to 
decrease along this second branch. Beyond B the mass increases again, and soon rises steeply along this third branch. 
But the radius (soon) starts to increase as well, and then rises further along this branch, almost up to the cusp C. 
Thus while the boson stars are already compact objects on the first branch, they become even much more compact 
objects along the third branch. 

This behavior is very reminiscent of the phases of compact fermionic stars, where the lower density stars represent 
white dwarfs, while the higher density stars represent neutron stars or their variants. The neutron stars seize to exist, 
when the black hole limit is reached. Close to this limit, the solutions exhibit a spiralling behavior as seen in the 
schematic drawing in Fig. [T] And indeed, when we include the black hole limit, given by the Schwarzschild relation 
2GM = R, in the figure for the boson stars, Fig. [6] (right), we observe that the spiral is precisely formed, when the 
family of boson star solutions approaches the black hole limit. 
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Figure 6: A comparison of the two definitions Eq. (|27[) for the boson star radius for the family of solutions with k — 0.002. The 
cusps A through D are also marked (left). The mass versus the boson star radius Ri for the families of solutions with k = 1.0, 
n = 0.2 and n — 0.002 together with the corresponding Schwarzschild black hole curve (right). 



Finally, in Fig. [7] we exhibit the mass of the boson stars versus the central particle number density J*(0) (left) and 
versus the central energy density p c (Q) (right). Clearly, the first branch has much lower central density than the third 
branch, just like the white dwarf branch and the neutron star branch of the compact fermionic stars [l| . 
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Figure 7: The mass versus the central particle number density j (0) (left) and the central energy density p c (0) (right) of the 
families of boson star solutions with k = 0.001 and k = 0.1. 



4- Stability analysis 



Let us now address the stability of these families of spherically symmetric boson stars. The above observed analogy 
between the branches of compact stars and the branches of boson stars, suggests to begin by recalling the stability 
properties of the compact stars as indicated in Fig. [TJ following Shapiro and Teukolsky |l|. 

The white dwarf branch ending at the cusp A is stable. It is followed by an unstable branch, ending at the cusp 
B. The neutron star branch from B to C is again a stable branch. The spiral beyond C is unstable. In particular, 
one mode becomes unstable at A, and turns stable again at B. At C the mode becomes unstable again, and at each 
following extremal point in the spiral another mode turns unstable. Thus there are two physically relevant stable 
branches, the lower density white dwarf branch and the higher density neutron star branch. The latter ends close to 
the black hole limit. 

For mini-boson stars a mode analysis has been performed by Lee and Pang [19j to determine the stability of the 
solutions. Mini-boson stars possess only a single stable branch followed by the spiral. As in the case of the compact 
fermionic stars, at each following extremal point in the spiral another mode of these boson stars turns unstable. 
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On the other hand, the stability of boson stars has been analyzed by envoking the arguments of catastrophe theory 
[20l l29l . [30} ) . We briefly recall the procedure employed by Tamaki and Sakai (29|, [3(| , applying it to the above family 
of boson star solutions. 

An essential point in utilizing catastrophe theory is to select an appropriate set of behavior variable(s) and control 
parameter (s). A behavior variable should be a quantity that describes the behavior of the system uniquely when the 
control parameters change their values. Following Tamaki and Sakai [H, we choose the coupling constant n and 
the charge Q as the two control parameters, and we choose the frequency lo s as the single behavior variable. (In 
the relevant range of solutions up to the cusp C, the variable uj s is indeed unique, but not beyond. In contrast, the 
variable cf>Q is unique in the spiral beyond the cusp C, but is not necessarily unique in the vicinity of C right before 
the cusp.) 

To analyze the stability of these boson stars, we start from the equilibrium space A4 = {uj s ,k, Q}, exhibited in 
Fig. [2] According to catastrophe theory, the stability changes only at the turning points, where -^j- = (for fixed 
values of k). Therefore 

• we determine the turning points, where stability changes, denoting them by A, B, C, . . . ; 

• we plot the values of the turning points A, B, C, . . . versus k to obtain the control space C = {k, Q}, to read 
off the regions of stability, respectively, instability. 

We list the turning points in Table HI] together with the extrapolated values of the limiting solutions Qu m at the 
centers of the spirals [171 ] . 



K 


Qa 


Qb 


Qc 


Qd 


Qe 


Qf 


Qlim 


0.0002 


1174.78 


40.2495 


3.117 • 10 7 


2.040 • 10 7 


2.220 • 10 7 


2.181 • 10 7 


2.181 • 10 7 


0.0020 


377.455 


40.0820 


864746 


553687 


605046 


593437 


5.954 ■ 10 5 


0.0100 


169.834 


39.9989 


60416.2 


36845.6 


40849.6 


39949.6 


4.011 ■ 10 4 


0.0200 


120.476 


38.7401 


17702.5 


10350.7 


11593.1 


11307.0 


1.136 • 10 4 


0.1000 


52.22 


33.91 


684.741 


331.677 


391.874 


378.274 


380.4600 


0.2000 


35.763 


28.668 


134.72 


61.19 


72.96 


70.463 


70.8330 


0.4000 


23.66 


21.47 


28.89 


15.85 


17.85 


17.46 


17.5180 


1.0000 


12.22 






5.83 


6.31 


6.29 


6.2114 



Table II: Values of the charge at the turning points A, B, C, . . . , and of limiting solution Qu m at the center of the spiral for 
the fundamental boson stars for several values of the coupling constant k. As indicated by the line, the cusps B and C have 
merged and disappeared between k = 0.4 and « = 1. 

Each point of the eqilibrium space A4 = {uj s ,k, Q} represents a boson star configuration. The set of turning 
points partitions this phase portrait into subareas. In catastrophe theory, "stability" means stability with respect 
to local perturbations. According to catastrophe theory, passing a turning point means changing the stability of the 
boson star configurations. The solutions between two turning points then form branches, where all the configurations 
possess the same kind of stability (or instability). Thus there are S- and [/-branches representing stable and unstable 
configurations [4lT |42|. 

This reveals the strength of catastrophe theory: stability changes exclusively when passing a turning point, while 
all configurations of the considered system between two turning points possess the same kind of stability (or insta- 
bility). Thus, it is sufficient to consider only a single configuration of a branch to know the stability of all the other 
configurations of the same branch with respect to local perturbations. 

For boson star solutions we see a cusp-catastrophe because (leaving the spiral aside) their equilibrium space M. 
shows the characteristics of a Whitney-surface [28|, l43f . Therefore we conclude that the branch from the vacuum 
solution to the first turning point Qa is a stable one |28l. 130} . The next branch from Qa to is unstable, and the 
branch from Qb to Qc is stable again 30]. This comprises a complete Whitney surface with a stable upper and 
lower sheet and an unstable area inbetween. The next turning points are part of the spiral. Such spirals are described 
by Arnold [271 ] as "limit cycles" , where the equilibrium states lose their stability. Then all configurations, which are 
part of the spirals, are unstable. We note, that this analysis is in complete agreement with the discussion above for 
compact fermionic stars. 

Let us finally exhibit the stability properties of the boson stars by presenting the control space C = {k, Q} in Fig. [5] 
The control space is the projection of the catastrophe map x{M) into the control plane. In the regions denoted by Si, 
SiU (i — 1,2) and N, there is one stable solution, there are i stable solution(s) and one or more unstable solution(s), 
and there are no equilibrium solution, respectively. The area delimited by the turning points Qc and Qd contains 
beside one stable solution all the unstable solutions of the spiral. 
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The control space reveals that there are two areas, where only a single stable solution exists. The first area is 
associated with the stable branch terminating at A. This is the branch formed by the stable lower density boson 
stars. The second area is associated with the branch from B to (almost) C. This branch comprises the stable high 
density boson stars. As the coupling constant k increases, the width of the second area diminishes and shrinks to zero 
at the critical value of n. This transition is consistent with the fact that for large n the solutions tend to mini-boson 
star solutions, which possess only a single stable branch. 




Figure 8: Control space C = {k, Q} for the fundamental boson stars. In the areas denoted by Si, SiU (i — 1, 2) and N, there is 
a single stable solution, there are i stable solution(s) and one or more unstable solution(s), and there is no equilibrium solution, 
respectively. The axes are semi-logarithmic (left) and double logarithmic (right). 



B. Rotating Boson Stars 



Let us now turn to rotating boson stars. These stationary axially symmetric configurations have a finite angular 
momentum proportional to their particle number, J = nQ (Eq. 122 j). with the "quantum number" b / 0. In the 
following we extend our above analysis of the physical properties of boson stars to the rotating case, focussing on 
solutions with n = 1 and positive parity, i.e. n = 1 + boson stars, that were obtained previously [17l Il8j. 



1. Equilibrium space of n = 1 boson stars 

As in the case of the spherically symmetric boson stars, we begin our analysis of the n = 1 + rotating boson stars 
by considering the equilibrium space A4 = {lu s , k, Q}. Again, for a fixed value of the dimensionless coupling constant 
k a family of stationary solutions exists in the frequency range ojo(k) < uj s < w max , where uiq(k) increases with n, 
whereas uj max is still given by the boson mass [TtJ ■ This is seen in Fig. where we illustrate the equilibrium space 
A4 = {uj s , k, Q} of the n = 1 + boson stars. 

Clearly, the structure of the equilibrium space Ai = {lo s , k, Q} of the n = 1 + boson stars is very similar to the one 
of the n — solutions. All families of solutions (for fixed k) start from the vacuum Q = at the upper limit of the 
frequency, and all form spirals at the lower end. However, we note that as n increases, the n = 1 + spirals become 
more elongated in this equilibrium space with respect to w s than their n = counterparts. 

We could also consider the set of equilibrium solutions with respect to other sets of parameters. For the n = boson 
stars we considered the alternative equilibrium space M — {4>q, M, k}. Here we would have to replace the variable (f>o 
by another variable, however, since for rotating boson stars 4>o — 0. Such an alternative variable for n = 1 + boson 
stars would be = <j>'(fi) 0]. 



2. Cusp structure of n = 1 boson stars 



Let us next turn to the cusp structure of these rotating boson stars. As in the case of the n = solutions, for not 
too large values of the coupling constant k a first branch of rotating boson star solutions resides between the vacuum 
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Figure 9: Equilibrium space M = {uj s , k, Q} for rotating boson stars with n — 1 + in the range 0.04 < k < 0.3. 

M = Q = and the local maximum of the mass and charge, Ma, Qa, at the first cusp A. Next, a second branch 
connects the first cusp A with the second cusp B, where the mass and charge have a local minimum, Mb, Qb- Then 
a third branch emerges from the second cusp B and extends up to the third cusp C, where the mass and charge reach 
their global maximum, Mq, Qc- Beyond C, finally, a spiral is formed. We list these turning points in Tab. IIIII 





Qa 


Qb 


Qc 


Qd 


0.04 


226.888 


135.570 


5121.20 


1630.19 


0.06 


184.8 


130.123 


2408.91 


434.076 


0.10 


140.631 


118.113 


910.699 


77.434 


0.20 


96.197 


94.1 


237.173 


35.984 


0.30 






110.958 


23.462 



Table III: Turning points for rotating boson stars with n = 1 for several values of the coupling constant k. 

Analogous to the non-rotating case, the values of the mass Ma and the charge Qa at the cusp A increase, ass^O. 
Indeed, we observe the same kT 1 / 2 dependence, as seen in Fig. [10] The values of the mass Mb and the charge Qb 
at the cusp B, on the other hand, are expected to tend to constant values in the limit k — > 0, corresponding now 
to the values of the mass and the charge of the unique minimum of the rotating Q-ball solutions in flat space-time. 
(The available data do not yet suffice to fully demonstrate this behaviour, however.) The values of the mass Mc 
and the charge Qc at the cusp C increase again without bound as k — > 0. The available data do not yet exhibit the 
expected limiting k~ 3 / 2 dependence, but a somewhat deviating n dependence. This is analogous to the non-rotating 
case in this range of values of the coupling constant n, where the limiting k -3 / 2 dependence is only approached for 
considerably smaller values of n. 

For larger values of n, however, we observe a somewhat different pattern for the cusps A, B and C in the case of 
the rotating boson stars as compared to the non-rotating case. In the rotating case, the extrema A and B merge, and 
disappear at a critical value of k, while C remains, as seen in Tab. IIIII and Fig. [TU] In the non-rotating case it is B 
and C , which merge, while A remains. In both cases, however, for large values of n the solutions exhibit the same 
pattern as the mini-boson stars: they possess a single physically relevant branch, beyond which a spiral is formed. 

3. Size of n — 1 + boson stars 

As in the case of the n = boson stars, we would like to study the compactness of the rotating boson stars and 
their proximity to the black hole limit. For that purpose we consider several possiblities to obtain a measure for the 
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Figure 10: Coupling constant n versus the mass M (left) and versus the particle number Q (right) for the cusps A, B and C, 
together with functions approximating the behavior of A and B in the limit k — >• 0. Note, that for the larger values of k the 
cusps A and B are no longer present. 



size of these axially symmetric boson star solutions. We define area type radii via introducing first the function R(r) 



1 



1/2 



R ( r ) = [ ^ /. ) , (28) 



and then using this R(r) to obtain as measures for the size 

Q 



Ra = [WbEii?. ( 29) 



and 



. 2 _!f{R{r)Y^—g<Pr 



Ra2 = J " q V • (30) 

Similarly, we define circumferential type radii via first introducing the function 

R ( r ) = hj c ^ (31) 

and then using this R(r) to obtain Rq as in Eq. (|29[) and R 2 C2 as m ^q. (|30| . All these definitions give rather similar 
results. We demonstrate this (in part) in Fig. [TTJ where we compare the radii Ra2, Ra an d Rc2- 

The figure clearly shows, that the dependence of the mass on the size is very similar for rotating and non-rotating 
boson stars. Focussing on the smaller values of k (away from the mini-boson star limit) the mass increases along 
the first branch, while the size decreases, until the cusp A is reached, marking the end of the lower density phase. 
Between A and B the mass decreases with decreasing size. Beyond B the mass then rises steeply again, while the 
size also increases along (most of) this third branch. This high density phase of the rotating boson stars then ends 
at the cusp C, where the spiral starts. 



4- Black hole limit 



As in the case of the non-rotating boson stars, we would now like to address the black hole limit for these rotating 
boson stars. For this purpose we employ these radii to compare the boson star masses and sizes with those of the 
corresponding rotating Kerr black holes. In particular, for the set of values of the size Ra2 and angular momentum 
J of the boson stars, their mass M is compared with the mass Mbh of the corresponding Kerr black holes, which 
are evaluated for the same values of Ra2 an d J , where Ra2 is defined via the event horizon area. We exhibit this 
comparison in Fig. [T^] for k = 0.04. We note, that because of this construction of the Kerr curve, it has as many 
branches as the boson star curve. Interestingly, we observe that as in the case n = 0, also the rotating boson stars 
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Figure 12: Comparison of the boson star mass M versus the size Ra2 for k = 0.04 with the mass Mbh of the Kerr black hole 
obtained for the same size R and angular momentum J (left); comparison of the ratio J/M 2 versus the ratio Ra2/M of the 
boson stars (k = 0.04) with the Kerr black hole curve (right). The point e indicates the onset of an ergoregion. 



are very close to the black hole limit, when the spiral is formed, as seen in Fig. [12] (left). (The spiral is enlarged in 
the inset, which also features the Kerr black hole branch corresponding to the steep high density boson star branch.) 

For a given mass Mbh, the angular momentum Jbh of a Kerr black hole cannot exceed a certain bound, the Kerr 
bound, Jbh /Mbh — 1- The Kerr bound is saturated for extremal black holes. Higher angular momenta correspond 
to naked singularities. When considering this scaled angular momentum Jbh/Mq H versus the scaled area Abh/M§ h 
of the black holes, all Kerr black holes fall onto a single line, starting from the point corresponding to the set of 
Schwarzschild solutions and ending at the point corresponding to the set of extremal Kerr solutions. The upper part 
of this line representing fast rotating Kerr black holes is exhibited in Fig. (right). 

In our comparison of the rotating boson stars with the Kerr black holes we would now like to know, how close the 
boson stars are to this Kerr bound, when they are highly compact. Therefore we also show J/M 2 versus Rai/M for 
such a set of highly compact rotating boson stars in Fig. Q2] (right). In particular, we have also indicated the point 
C of the boson star curve, where the spiral starts, and the point D, which is located inside the spiral. We note that 
the family of boson stars assumes at C its minimum value of J/M 2 . For k = 0.04, the case shown in the figure, 
J/M 2 k, 0.7 at C. When k increases, the value of J/M 2 at C increases (k = 0.06: 0.78, k = 0.1: 0.86, k = 0.2: 0.92, 
k = 0.3: 0.94). Thus these highly compact boson stars close to the black hole limit can rotate with angular momenta 
close to the Kerr bound. 

In Fig. Q2] (right) we have also indicated the point, beyond which the rotating boson stars possess ergoregions 
[l8j . For k — 0.04, the formation of ergoregions arises only inside the spiral. For larger values of k the ergoregion 
formation starts already on the high density branch shortly before C is reached. We therefore conclude that for 
these n = 1 + rotating boson stars there appears to be a correlation between approaching the black hole limit and 
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developing an ergoregion. 



5. Stability analysis of n — 1 boson stars 



To address the stability of these rotating boson stars, we turn to the control space C = {k, Q}, presented in Fig. [T3] 
Since both n = and n = 1 + boson stars have analogous equilibrium spaces M. = {uj s ,k, Q}, the corresponding 
discussion of their stability is also analogous. From the analysis via catastrophe theory we conclude that there are two 
areas which possess only a single stable solution. The first area is associated with the stable branch terminating at 
A, formed by the stable lower density boson stars. The second area is associated with the branch from B to (almost) 
C, comprising the stable high density boson stars. As the coupling constant k increases, the width of the second area 
diminishes and shrinks to zero at a critical value of K. The width of the area labeled to represent 2 stable solutions 
as well as unstable solutions also diminishes with increasing k. The width shrinks to zero at the critical value of k, 
where the turning points A and B merge and disappear, while only C remains (apart from the critical points of the 
spiral). 
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Figure 13: Control space C = {«, Q} for rotating boson stars with n = 1 + in the range 0.04 < k < 0.3. 



According to catastrophe theory arguments the branch from the vacuum to the first turning point Qa is a stable 
one. The next branch from Qa to Qb is unstable, while the branch from Qb to Qc should be stable again. However, 
we must also take into account the emergence of ergoregions, which are associated with their own kind of instability 
[451 - I481 ] . While for small values of k ergoregions arise only beyond C somewhere in the spiral, which is unstable 
anyway, for larger values of k (k > 0.1) ergoregions appear already along the higher density branch from B to C, 
signalling an instability of these highly compact rotating boson stars close to C [3, |48| . 

Apart from this new type of instability, present only in the rotating case, the boson star solutions for the non- rotating 
(ri = 0) and the rotating (n = 1 + ) case exhibit a similar general pattern. 



IV. CONCLUSIONS 



We have addressed the physical properties of non-rotating and rotating boson stars, obtained with a self-interaction 
potential of the scalar field, which allows for non-topological soliton solutions in the absence of gravity. Such a self- 
interaction potential is crucial to find a much richer set of solutions than the ones obtained with only a mass term 
(mini-boson stars) and with a repulsive |$| 4 interaction. 

In particular, we note that there are two stable regions in the equilibrium space of the boson star solutions, when 
this solitonic self-interaction potential is employed, whereas there is only one such stable region for boson stars which 
do not possess a flat space-time limit. In fact, there is an interesting analogy with compact stars, which possess a 
lower density phase, the white dwarf phase, and a high density phase, the neutron star phase, since the boson stars 
also exhibit a lower density phase and a high density phase. Moreover, beyond the neutron star phase the compact 
stars exhibit an unstable spiralling phase very close to the black hole limit. Such an unstable spiralling phase very 
close to the black hole limit is also seen for the boson stars beyond their stable high density phase. 

A stable mini-boson star is an equilibrium state, where the Heisenberg uncertainty principle ArAp ~ irk provides 
the means to balance gravity and avoid collaps, below a (small) critical mass of the mini-boson stars [16j . In boson 
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stars with a repulsive |<I>| 4 potential term this self- interaction allows for much larger stable boson stars [10(. The 
solitonic self-interaction potential, on the other hand, has repulsive and attractive components, which dominate the 
features of the solutions in a large region of the equilibrium space, leaving here for gravity only a minor role to play. 

Indeed, for not too large values of the coupling constant k, the properties of the solutions follow rather closely those 
of the corresponding non-topological solitons. Only at the boundaries of the domain of existence, gravity becomes 
dominating. Here in flat space-time the solutions would grow without limit. Thus the single infinitely long stable 
branch of the non-topological solitons is reflected in the finite (mostly) stable branch B — C of the boson stars. Along 
(most of) this branch the radius increases as the mass increases. This branch ends, when the boson star's compactness 
approaches the black hole limit. The soliton solutions would simply cross this limit. But as long as the coupling to 
gravity is finite, no matter how small it is, this limit cannot be exceeded. Thus collapse is unavoidable, as signalled 
by the formation of a spiral in the equilibrium space of these stationary solutions. 

The single infinitely long unstable branch of the non-topological solitons, on the other hand, is reflected in the finite 
unstable branch A — B of the boson stars. But gravity allows for an additional stable branch — A, present even 
in the case of no self-interaction. Here the size of the solutions decreases as the mass increases, which reveals the 
dominating influence of gravity for this branch. 

When comparing the non-rotating (n = 0) and the rotating (n = 1 + ) boson stars we note, that many of their 
features are very similar. While the high density soliton-type boson star branch is bounded by the Schwarzschild 
black holes in the n — case, it is bounded by the Kerr black holes in the n = 1 + case. Because the centrifugal force 
will counteract the gravitational force, it is expected that rotation stabilizes a boson star [l6| , and the available data 
indeed show, that in the presence of rotation for given values of k higher masses are reached for the rotating boson 
stars. However, rotation also comes with a caveat for stability, since for globally regular objects such as boson stars 
the presence of an ergoregion implies an instability, associated with superradiant scattering [45T - |48j . 

Thus rotating boson stars become unstable, when they develop an ergoregion. For the n — 1 + boson stars considered, 
the ergoregion is formed either inside the spiral or very close to the cusp C. Therefore the possible presence of 



ergoregions put forward by Cardoso et al. 48( as an argument to exclude boson stars and various other black hole 



doubles as potential horizonless candidates for compact dark astrophysical objects, is not yet compelling. It remains 
to be seen, whether boson stars with appropriate values of the physical parameters to fit observational data will or 
will not suffer from such an ergoregion instability. Our analysis shows, that there are stable highly compact boson 
stars close to the black hole limit, that are at the same time close to the Kerr bound J/M 2 = 1. 

Finally, we would like to address the case of rotating boson stars with higher rotation quantum numbers n > 2. 



Such boson stars have been addressed by Ryan [13[ . Since the numerical analysis for these systems is difficult [17|, \li 



we have not yet accumulated a sufficient set of solutions to give a similar discussion of their properties as for the n = 
and n = 1 + boson stars. Still, we present in Fig. [14] a part of the equilibrium space A4 = {uj s , k, Q} of the n = 2 + 
boson stars. Whereas still higher values of the mass are reached on the soliton-type branch, we note, that the onset 
of the ergoregion instability happens earlier, thus decreasing the viablity range of these solutions. 




Figure 14: Equilibrium space M = {cj s ,k, Q} for rotating boson stars with n = 2 + in the range 0.04 < k < 0.3. The onset of 
the ergoregion is indicated by the asterisks. 



Clearly, further calculations of fast rotating boson stars are called for, not only in the solitonic model but also in 
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|$| -type models, where a possible ergoregion instability was not yet considered. Most interesting will, however, be 
an extension of the investigation of the dynamical evolution of boson stars [4l|, |42[ to the rotating case. 
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Appendix A: Construction of the solutions 



1. Boundary conditions 



The choice of appropriate boundary conditions must guarantee that the boson star solutions are globally regular 
and asymptotically flat, and that they possesses finite energy and finite energy density. 

For spherically symmetric boson stars boundary conditions must be specified for the metric functions f(r) and Z(r) 
and the scalar field function <f>(r) at the origin and at infinity. At the origin one finds the boundary conditions 

<9r/|r=0 = , d r l\ r = = , d r (f>\ r =o = . (Al) 

Note, that for spherically symmetric boson stars the scalar field has a finite value 4>o at the origin, 

cj>{r) = 0o + 0(r 2 ) . (A2) 

For r — > oo the metric approaches the Minkowski metric rj a ^ and the scalar field assumes its vacuum value $ = 0. 
Accordingly, we impose at infinity the boundary conditions 

f\r— >oo — 1 7 l\i — >oo — 1 ; 0| r— >oo = . (A3) 

For rotating axially symmetric boson stars appropriate boundary conditions must be specified for the metric func- 
tions f(r,0), l(r,9), h(r,6) and oj(r, 6), and the scalar field function <p(r,9) at the origin, at infinity, on the positive 
z-axis (9 = 0), and, exploiting the reflection symmetry w.r.t. 9 — > it — 9, in the xy-plane (9 = it/2). At the origin we 
require 

<3 r /| r=0 = , d r l\ r=0 = , h\ r=0 = 1 , uj\ r=0 = , </>| r=0 = . (A4) 
At infinity the boundary conditions are 

f\i — s-oo = 1 i l\r— yoo = 1 ) h\ r — >oc = 1 ? w |t — s-oo — , 0|r— ¥00 — , (A5) 

and for 9 = and 9 = tt/2, respectively, we require the boundary conditions 

d s f\e=o = , d e l\ g=0 = , h\ g=0 = 1 , dguj\ e = = , 0| e=o = , (A6) 
and for even parity solutions 

def\e=Tr/2 = , d e l\ g=7: /2 = , d e h\e^/ 2 = , dguj\ e ^/ 2 = , <9e0| e=7r / 2 = , (A7) 
while for odd parity solutions <p\e=Tx/2 = 0. 



2. Numerical methods 



First of all, because of the power law fall-off of the metric functions, we compactify space by introducing the 
compactified radial coordinate 

(A8) 



l + r 

Then the resulting set of equations is solved numerically subject to the above boundary conditions. 



17 



For spherically symmetric non-rotating solutions (n = 0) the set of equations depends only on the radial coordinate. 
It is solved numerically by employing a collocation method for boundary-value ordinary differential equations devel- 
oped by Ascher, Christiansen and Russell Here the damped Newton method of quasi- linearization is applied. At 
each iteration step a linearized problem is solved by using a spline collocation at Gaussian points. 

Rotating solutions are obtained when n ^ 0. The resulting set of coupled non- linear partial differential equations is 
solved numerically by employing a finite difference solver [50j, based on the Newton-Raphson method. The equations 
are discretized on a non-equidistant grid in f and 9. Typical grids used have sizes 90 x 70, covering the integration 
region 0<f<land0<6»< tt/2. 
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